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Abstract 

We report the measurement of the direct CP asymmetry in the radiative B —» Xs+d') decay using a data sample of 
(772 ± 11) x 10 6 BB pairs collected at the Y(4S) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ 
collider. The CP asymmetry is measured as a function of the photon energy threshold. For E] >2.1 GeV, where E*, is the 
photon energy in the center-of-mass frame, we obtain Acp{B — >• A's+dy) = (2.2 ± 3.9 dh 0.9)%, consistent with the Standard 
Model prediction. 

PACS numbers: 11.30.Er, 13.25.Hw 


The radiative electroweak transitions b —> s'y and 
b —> d'y proceed via flavor-changing neutral currents in¬ 
volving loop diagrams. These decays are sensitive to pos¬ 
sible contributions from new heavy particles occurring in 
the loop, which modify the branching fractions and CP- 
violating effects predicted in the Standard Model (SM). 
The decay rates, including QCD corrections, can be ex¬ 
pressed by an effective Hamiltonian and calculated using 
the Operator Product Expansion approach. In the lead¬ 
ing and next-to-leading order logarithmic approximation, 
the branching fractions and CP asymmetries are pro¬ 
portional to the dipole operators P 7 and [l|. New 
physics effects would modify the corresponding Wilson 
coefficients C 7 and C&. 

The CP asymmetry ( Acp ) in B —> Xs+d'y decays is 
defined as: 


, ,5 „ T(B ->• X s+d7 ) - T(B -)• X- a+ n) 

cp{ s+dl) ~ r {B X s+dl ) + T(B -+ X s+ry )’ 

(1) 

where T(B — > X s+d "f ) represents the decay rate of the 
B° or B~ meson into the radiative final state. In the 
following, charge-conjugate states are included implic¬ 
itly. The X s+ d states represent all possible hadronic 
final states derived from b —> sj or b —> dj transi¬ 
tions. The SM predicts Acp for the these two transitions 
in the ranges —0.6% < Acp(B -A X s y) < 2.8% and 
—62% < Acp(B —>■ X d j) < 14% Q. Even though the 
individual CP-violating effects could be large, the CP- 
violating contributions cancel when both are considered 
inclusively due to CKM unitarity, and the theory errors 
cancel almost perfectly except for small [/-spin breaking 
corrections [§i], additionally, the inclusive asymmetry is 
insensitive to the choice of photon energy cutoff Q. This 
precise SM prediction of Acp (B —> X s+ d 7 )= 0 , serves 
as a clean test for new CP-violating phases acting in the 
decays. New physics (NP) scenarios such as supersym¬ 
metric models with minimal flavor violation predict Acp 
(.B —>• Xg+d'y') up to a level of +2 %. In more generic NP 
scenarios, the asymmetries Acp (B —> X., 7 ) and Acp 
(.B —>• Xd'y) do not cancel and Acp (B —> X s+d 'y) is the 
most sensitive observable, with values as large as 10% Q. 

Previous measurements of Acp (B —> X s+ d 7 ) have 
been performed by CLEO @| and BaBar [g| and are sta¬ 
tistically limited. Belle has performed a measurement 
of the inclusive branching fraction .7J. The asymmetry 
Acp (B —> X s "f) has been measured separately as the 
sum of exclusive decays @, 0 ]. In this letter, we present 


the first Belle measurement of Acp(B — > X s+ d"f). We 
profit from the large data sample available at Belle to 
achieve a higher statistical precision. 

The states X s+d include resonant contributions such 
as K*( 892), p and to, and non-resonant contributions. In 
order to be sensitive to all X s+d states, the selection is 
based on the high-energy-photon signature of the transi¬ 
tion, i.e. the radiated photon is the only reconstructed 
particle from the B —> X s+d j decay. While this ap¬ 
proach does not exclude explicitly possible contributions 
from B —>• Xc'y or B —> X u c decays, such contributions 
are very small in the SM 10] and will be neglected in 
this analysis. To tag the signal B flavor, we use the fact 
that B mesons are produced in pairs from the reaction 
e + e _ —> Y(4 S) —> BB. The flavor of the signal B me¬ 
son is determined by tagging the flavor of the other B in 
the event, using a charged lepton (e, p) consistent with 
the semileptonic decay of the other B. The B flavor and 
lepton charge in semileptonic decays are directly related. 

Since the expected CP violation is very small and pre¬ 
cisely calculable, all effects that could bias the measure¬ 
ment must be carefully quantified. A measurement bias 
is introduced if the selection procedure, track reconstruc¬ 
tion, or particle identification favors a particular charge. 
These effects are quantified in different control samples. 
In this analysis, we also test the independence of Acp 
with respect to the choice of cutoff energy, by measuring 
it as a function of the photon energy threshold. 

This analysis uses the 711 fb _1 sample recorded at the 
T(4S) resonance by the Belle experiment at the KEKB 
storage ring jTlj, containing (772 ± 11) x 10 6 BB pairs. 
An 89 fb -1 sample recorded at a center-of-mass (CM) en¬ 
ergy 60 MeV below the resonance is used to study con¬ 
tinuum background (e + e~ —> qq 1 where q = u,d,s,c ); 
the former sample is denoted on-resonance and the lat¬ 
ter off-resonance. The Belle detector is a large-solid- 
angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cerenkov counters 
(ACC), a barrel-like arrangement of time-of-flight scintil¬ 
lation counters (TOF), and an electromagnetic calorime¬ 
ter comprised of CsI(Tl) crystals (ECL) located inside 
a super-conducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside of the 
coil is instrumented to detect mesons and to identify 
muons (KLM). The detector is described in detail else¬ 
where [r|j]. 
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Monte Carlo (MC) simulation samples were gener¬ 
ated to study continuum background, BB decays and 
B — > Xg'y signal events. The size of the BB MC sample 
is equivalent to ten times the integrated luminosity of the 
data. The size of the continuum MC sample corresponds 
to the integrated luminosity of the on-resonance sample. 
The generation of the signal B —>• X s 7 decays follows the 
theoretical prediction of the Kagan-Neubert model [T3| 
with parameters mb = 4.574 GeV and p% = 0.459 GeV 2 
representing the 6 -quark mass and mean kinetic energy. 
The signal sample contains 2.6 million B —» X s 7 events, 
which corresponds roughly to five times the number ex¬ 
pected in data. The BB and B — > X s y MC samples 
included B°-B° mixing. 

In this analysis, tracks passing very far from the inter¬ 
action point or compatible with a low-momentum par¬ 
ticle reconstructed multiple times as it spirals in the 
CDC are excluded. For photons, minimum energies of 
100 MeV, 150 MeV and 50 MeV,respectively, are re¬ 
quired in the forward, backward and barrel regions of 
ECL, defined in Ref. [12| . These requirements suppress 
low-energy photons resulting from particle interactions 
with detector material or the beam pipe. All particles 
are used to calculate kinematic and topological variables. 

The signal photon candidates are selected as connected 
clusters of ECL crystals in the polar angle 32.2° < 9 1 < 
128.7° with a CM energy 1.4 GeV < E* < 4.0 GeV. The 
polar angle is measured from the z axis that is collinear 
with the positron beam. The ratio of the energy deposit 
in the central 3x3 crystals to that in the central 5x5 crys¬ 
tals must be larger than 90%. Photons from the decays 
7r° ( 77 ) —> 77 are rejected by using a veto based on the pho¬ 
ton energy, polar angle and the reconstructed diphoton 
mass, as described in Ref. 0 - The signal region includes 
photons with CM energy 1.7 GeV < E* < 2.8 GeV; the 
sidebands E* < 1.7 GeV and E* > 2.8 GeV are used 
to study the normalization of BB and continuum back¬ 
ground components, respectively. 

The lepton candidates used for tagging (tag lepton) are 
reconstructed as tracks in the SVD and CDC. We limit 
the impact parameters along the z axis to \dz\ < 2 cm 
and dr < 0.5 cm, require at least one hit in the SVD, 
and choose a momentum range in the CM frame of 
1.10 GeV < p} < 2.25 GeV. The upper-momentum 
bound reduces continuum background as it is near the 
kinematic limit for leptons from B decays. The lower 
bound ensures that most of the selected leptons origi¬ 
nate directly from a B meson, which is important for 
flavor tagging. Electron candidates are identified by con¬ 
structing a likelihood ratio based on the matching of the 
cluster in the ECL and the extrapolated track, the ratio 
between its energy and momentum, the shower shape in 
the ECL, the energy loss in the CDC, and the light yield 
in the ACC. The polar angle requirement for electrons is 
18° < 0 e < 150°. Muon identification uses a likelihood 
ratio determined from the range of the track and the nor¬ 


malized transverse deviations between the track and the 
KLM hits associated to it. The polar angle requirement 
for muons is 25° < 9^ < 145°. 

After this initial selection, the sample is dominated 
by continuum background, which constitutes 77% of the 
total yield; the signal component amounts only to 1 % 
as can be seen in Fig. QJa). To suppress the contin¬ 
uum background, we use a Boosted Decision Tree (BDT), 
that is trained to achieve the best discrimination between 
continuum and signal events. Eighteen kinematic, event 
shape and isolation variables are used as input for the 
BDT: eleven Fox-Wolfram moments 0 , constructed in 
three sets in which ( 1 ) all particles in the event are used, 
(2) the signal photon is excluded and (3) both signal pho¬ 
ton and tag lepton are excluded; the magnitude and di¬ 
rection of the event’s thrust vector; the distance between 
the photon cluster and the closest extrapolated position 
of a charged particle at the ECL surface; the angle be¬ 
tween the directions of the photon and tag lepton; the 
RMS width of the photon cluster; the scalar sum of the 
transverse momenta of all reconstructed particles; and 
the square of the missing four-momentum, calculated as 
the difference between the total beam energy and the mo¬ 
menta of all reconstructed particles. The BDT is trained 
using continuum and B —>• X s ^ MC samples. The se¬ 
lection criterion on the BDT output classifier variable is 
chosen to minimize the expected statistical uncertainty 
on Acp ■ The BDT classifier distribution and selection 
criterion are shown in Fig. [2] The photon spectrum after 
continuum suppression is shown in Fig. [ljb) for MC and 
on-resonance data, in this plot we include statistical un¬ 
certainties and systematic uncertainties that come from 
calibration and normalization factors, that cancel in the 
measurement of Acp- 

After the selection, in the region 1.7 GeV < E* < 
2.8 GeV, we find 21400 (21608) events tagged with a 
positive (negative) lepton in the on-resonance sample and 
2623 ± 140 (2728 ± 143) events tagged with a positive 
(negative) lepton in the off-resonance sample. The off- 
resonance events are corrected as they have, on average, 
lower particle energies and multiplicities due to the lower 
CM energy. Additionally, the off-resonance yield is scaled 
to take into account the difference in luminosities and 
cross-sections. 

The signal fraction is 21.2% while the continuum back¬ 
ground fraction is 12.4%. The BB background contains 
photons from several processes. The dominant sources 
are photons from ir° —>• 77 decays, which make up 49.5% 
of the total yield and photons from 77 —» 77 , contributing 
7.9%. Photons from beam background are 2.2% of the 
total contribution. Electrons and hadrons misidentified 
as photons are small contributions of 0 . 8 % and 0 . 2 %, 
respectively. Other photons, mainly from decays of w, 
t)' and J/ip mesons, and bremsstrahlung, including final 
state radiation [l 6 |, comprise the remaining 5.8%. The 
B —>• V s+ d 7 signal is obtained by subtracting the contin- 
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FIG. 1: Photon energy spectrum in the CM frame showing on-resonance data, off-resonance data for continuum, and MC 
simulation. The spectrum is shown (a) before and (b) after continuum suppression. In (a), the MC signal is additionally plotted 
scaled a factor of fifty to show its expected position. In (b), the MC error includes statistical and systematic uncertainties 
coming from calibration and normalization factors that cancel in the measurement of Acp- 



FIG. 2: Output of the BDT. The continuum distribution cor¬ 
responds to off-resonance data. The vertical line denotes the 
minimum requirement on this variable. 


uum and BB contributions. The BB background com¬ 
ponents are calibrated using data, as described below. 
All corrections and calibrations applied to MC and off- 
resonance data are determined and performed indepen¬ 
dently of the tag charge. The subtraction of background 
is done for each charge individually. 

The rejection of events containing 7 r° or 77 will fail in 
cases where the decay is very asymmetric and the second 
photon has an energy below the threshold, making the 
reconstruction of the 7 r° or 77 impossible. To properly 
normalize these components, the veto is removed and, 
for each combination of the prompt photon with another 
photon in the event, the diphoton mass m 77 is calculated. 
A fit to the 7 T° and 77 masses is performed to estimate 
the number of these mesons in data and MC. The fit is 
performed in eleven meson momentum bins between 1.4 
and 2.6 GeV and the ratio of data to MC yields is used 


as a correction factor. 

Some background components have a non-vanishing di¬ 
rect CP asymmetry that could impact our measurement. 
Most have negligible contributions to the decay rate ex¬ 
cept for B —>• X s rj decays, which comprises 1.2% of the 
rate according to the MC prediction, with a branching 
fraction B{B X s r]) = (26.1 ± 3.0 ±3 ? x 10 " 5 and 
a CP asymmetry Agp(B —> X s rj) = (—13 ± 5)% mea¬ 
sured by Belle 17]. The MC is corrected to model this 
effect properly. 

The B —> X s+c i / y photon energy spectrum for positive 
and negative tagged events after subtracting all the back¬ 
ground is shown in Fig. [3] The measured asymmetry, 
_4™p as , is calculated using Eq. dTJ expressed in terms of 

N+ - N~ 

the charge-flavor correlation: = ——-. Here, 

6 CF N+ + N~ 

N + and N~ represent the total number of events tagged 
by a positive or negative lepton for a given photon en¬ 
ergy threshold. The energy thresholds range from 1.7 to 
2.2 GeV. 

The measured values must be corrected due to pos¬ 
sible asymmetries in the BB background that is sub¬ 
tracted (Mbkg) and possible asymmetries in the detection 
of leptons Adet • An additional correction arises from the 
probability that the reconstructed lepton has a wrong 
charge-flavor correlation, the so-called wrong-tag proba¬ 
bility (w). The corrected asymmetry is given by: 

Acp = y^(^cp s - Ak g - Aiet). (2) 

The correction Mdet accounts for a possible asymmetry in 
the identification efficiency between positive and negative 
charged leptons (Mlid) and a possible asymmetry be¬ 
tween the reconstruction of positive and negative tracks 
(Atrack)- Alid is determined using aB-> XJ/i/j (£ + £~) 
sample, where the selection efficiencies of positively and 
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FIG. 3: The photon energy spectrum in the CM system af¬ 
ter subtracting all the background, with vertical dashed lines 
showing the signal region. The positive tagged events are 
shown as circles and the negative as squares. Statistical and 
systematic uncertainties are included. 


negatively charged electrons and muons are determined 
by performing fits to the invariant dilepton mass mu 
for singly- and doubly-identified lepton candidates. The 


asymmetry is calculated as: Aid = 


This 


measurement is performed in the full kinematic region, 
in eleven laboratory-frame momentum bins and eight 
polar-angle bins. The asymmetries for electrons and 
muons are measured to be Aro(e) = (0.26 ± 0.14)% 
and Aid(aO = (—0.03 ± 0.03)%, and average to Aid = 
( 0.11 ± 0.07)%. The asymmetry Arack is measured with 
partially and fully reconstructed D* with D* —>• ttD°, 
D° —» 7T7T Kg, Kg —» 7r + 7r - decays, to be Arack = 
(—0.01 ± 0.21)%. The total detector-related asymmetry 
is Aet = ( 0.10 ± 0 . 22 )%. 

We measure Akg in the low-energy sideband E* < 
1.7 GeV. The asymmetries measured in data and MC 
are Akg(data) = (—0.14 ± 0.78)% and Ak g (MC) = 
(—0.26 ± 0 . 21 )%, which are consistent with zero within 
uncertainties. The asymmetry in the BB data is taken as 
a correction to Acp- Since this is an asymmetry in the 
BB background, the correction is proportional to the ra¬ 
tio of BB to signal events in the signal region, the ratios 
are taken from MC simulation. 


The wrong-tag probability has contributions from 
B°B° oscillations (uw), secondary leptons (w se c) and 
misidentihed hadrons (wmisro) and is given by u> = 
uJose + Wsec + WmisiD- The oscillation term is equal to 
the product of the mixing probability in the B°B° sys¬ 
tem Xd = 0.1875 ± 0.0020 [lj|, the fraction of neutral B 
mesons from the T(4S) decay, /oo = 0.487 ± 0.006 [l8j . 
and the fraction of leptons coming directly from a B de¬ 
cay, which is estimated to be 91.1% from MC, result¬ 
ing in w osc = 0.0832 ± 0.0015. Secondary leptons are 
true leptons that do not come directly from a B me¬ 


son but rather from one of its decay daughters. We 
find w sec = 0.0431 ± 0.0036; this value is estimated 
from MC and the error based on the precision with 
which the B — > DX and D —)• Xlv branching fractions 
are measured. Misidentihed hadrons give the smallest 
contribution and consist of 7 r and K mesons faking a 
muon and, to a lesser extent, an electron. The corre¬ 
sponding wrong-tag probability is estimated from MC, 
where the fraction of misidentihed hadrons is determined 
by studying D* + — > D 0 (K~n + )ir + decays. After ap¬ 
plying the same selection criteria for 7 r and K candi¬ 
dates as for tag leptons, the fraction of hadrons pass¬ 
ing the selection in the MC is corrected and we obtain 
WmisiD = 0.0069 ± 0.0034. The total wrong-tag probabil¬ 
ity value is w = 0.1332 ± 0.0052. 

The asymmetries Aet and Akg are the dominant 
uncertainties on Acp and are additive. An additional 
multiplicative systematic uncertainty arises from the 
wrong-tag probability, leading to a relative uncertainty 
A Acp/Acp = 0 . 01 , much less than the additive uncer¬ 
tainties. 

Finally, as some background events remain in the low 
energy range after subtraction, we scale the BB compo¬ 
nent to match the data yield below 1.7 GeV and recal¬ 
culate Acp • The difference between this value and the 
nominal is taken as an additional systematic uncertainty. 

In Table HI the measured and corrected values of Acp 
are summarized for 0.1 GeV steps in the E* threshold 
from 1.7 to 2.2 GeV, with a E* upper bound of 2.8 GeV. 
The statistical and systematic uncertainties are summa¬ 
rized in Table HT1 The statistical precision is improved in 
comparison to previous measurements it is, how¬ 

ever, the limiting factor in the measurement and is af¬ 
fected by the size of the continuum sample. As an exam¬ 
ple, for the 1.7 GeV threshold, the total 4.4% statistical 
uncertainty incorporates a 3.0% contribution from T(45') 
data and 3.1% from off-resonance data. The dominant 
systematic uncertainty arises from the asymmetry in the 
BB background. The asymmetry is consistent with zero 
across the different photon energy thresholds. 

In conclusion, we have measured the direct CP asym¬ 
metry Acp(B Xg+d'y)- The measurement is per¬ 
formed using (772± 11) x 10 6 BB pairs for photon energy 
thresholds between 1.7 and 2.2 GeV. As a nominal re¬ 
sult we choose the 2.1 GeV threshold since it has a low 
uncertainty and keeps a large fraction of signal events: 
Acp(B — > X s+dl ) = (2.2 ± 3.9 ± 0.9)%, consistent with 
the SM prediction. This is the first Belle measurement 
of this asymmetry and the most precise to date. 
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the Nil, and PNNL/EMSL for valuable computing and 
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from MEXT, JSPS and Nagoya’s TLPRC (Japan); ARC 
(Australia); FWF (Austria); NSFC (China); MSMT 
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TABLE I: CP asymmetry, in percent, for different photon energy thresholds, the E* kinematic limit is 2.8 GeV. For the 
measured asymmetry, only the statistical uncertainty is shown; for the corrected asymmetry Acp, statistical and systematic 
uncertainties are given. The ratio B/S representes the ratio in the number of BB to signal events that is used to scale the 
asymmetry _4bkg- An additional systematic uncertainty related to the wrong-tag probability is not explicitly listed but is taken 
into account in the total uncertainty; its relative value is 1%. The systematic contributions are added in quadrature. 


E* (thresh.) 

imeas 
^CP 

B/S 

*^bkg -^det 

MC stats. 

BB norm. 

Acp 

1.7 GeV 

1.3 ±3.1 

3.20 

-0.4 ±2.5 0.1 ±0.2 

±0.8 

±0.5 

2.2 ±4.3 ±3.5 

1.8 GeV 

2.0 ±3.0 

2.41 

-0.3 ±1.9 0.1 ±0.2 

±0.7 

±0.1 

3.0 ±4.1 ±2.7 

1.9 GeV 

0.9 ±2.9 

1.70 

-0.2 ±1.3 0.1 ±0.2 

±0.6 

±0.3 

1.4 ± 4.0 ± 1.9 

2.0 GeV 

1.6 ±2.8 

1.10 

-0.2 ±0.9 0.1 ±0.2 

±0.5 

±0.0 

2.2 ±3.8 ±1.3 

2.1 GeV 

1.6 ±2.9 

0.65 

-0.1 ±0.5 0.1 ±0.2 

±0.4 

±0.1 

2.2 ±3.9 ±0.9 

2.2 GeV 

1.1 ±2.9 

0.38 

-0.1 ±0.3 0.1 ±0.2 

±0.3 

±0.2 

1.4 ±3.9 ±0.6 


TABLE II: Absolute uncertainties in Acp, in percent. The systematic uncertainties are added in quadrature to yield the total. 


E* thresh. 

Statistical Total systematic Mdet Abkg 

MC stat. 

BB norm. 

Wrong tag 

1.70 GeV 

4.26 

3.52 

0.30 

3.40 

0.76 

0.42 

0.02 

1.80 GeV 

4.13 

2.72 

0.30 

2.56 

0.68 

0.53 

0.05 

1.90 GeV 

3.96 

1.92 

0.30 

1.81 

0.58 

0.10 

0.02 

2.00 GeV 

3.84 

1.32 

0.30 

1.17 

0.48 

0.19 

0.04 

2.10 GeV 

3.91 

0.86 

0.30 

0.70 

0.39 

0.12 

0.04 

2.20 GeV 

3.89 

0.59 

0.30 

0.41 

0.30 

0.04 

0.03 
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